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Method of creating images in computer tomography, and 
CT device for carrying out the method 

The invention relates to a method of creating images in 
5 computer tomography, in which an object under 
examination is scanned with the aid of a beam 
originating from a focus and a flat detector array, the 
output data determined is filtered in a suitable way 
and backprojected in order to obtain at least one slice 
10 which represents the absorption values of the section 
of the object under examination. Furthermore, the 
invention relates to a CT device which is suitable for 
carrying out this method. 

15 Methods of this type have become known under the term 
tt filtered backpro j ection" but, as 3D methods, that is 
to say in conjunction with a matrix-like detector 
array, do not furnish an image quality which is judged 
to be adequate in practice, since what are known as 

20 "cone beam artefacts" occur because of the conical x- 
ray beam. 

In addition, the disadvantage with this method is that 
redundant data, such as are produced during spiral 
25 scanning with low table advance as a result of multiple 
irradiation of one and the same voxel, is not used, 
which results in the radiation dose administered to the 
object under examination being used only incompletely 
for imaging. 

30 

Furthermore, there are thoughts relating to proceeding, 
in conjunction with 2D methods for image 
reconstruction, in such a way that preliminary images 
are calculated in large numbers by means of "filtered 
3 5 backpro j ection" from initial data which originates from 
sections of the focal path which are intrinsically 
inadequate for image reconstruction, the preliminary 
images being reformatted to form a final slice only in 
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a second step. These 2D methods are less useful for 
detector arrays with a large width, that is to say a 
great extent in the direction of the system axis, since 
an extremely large number of preliminary images then 
5 has to be processed, which is a problem even when high 
computing power is available. 

In all the aforementioned methods, the problem arises 
that, because of possible data redundancy which occurs 
10 during the scanning of the object under examination, 
image artefacts are produced which have a detrimental 
influence on the image quality. 

The invention is based on the object of specifying a 
15 method of the type mentioned in the introduction which 
permits the image quality to be increased. 

According to the invention, this object is achieved by 
the features of the independent patent claims. Special 

2 0 configurations and improvements emerge from the 

features of the subclaims . 

The inventors have discovered that in image processing 
it is firstly necessary to take account of the 

25 redundancies which occur with respect to a respectively 
considered voxel in the region under examination, in 
particular during the backpro j ection of the filtered 
data. Secondly, however, it is also necessary to take 
into account how the ray is positioned in the beam. For 

30 example, it has a decisively positive effect on the 
image quality if rays which are located centrally, 
based on the extent of the beam in the direction of the 
axis of rotation, are taken more into account when 
creating the image than rays which are positioned only 

3 5 marginally in relation to the beam, based on the same 

extent . 
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Accordingly, the inventors propose an improved method 
for creating images in computer tomography which has 
the following method steps: 

- to scan an object under examination with a beam 
5 originating from at least one focus and with a flat 

detector array having a large number of distributed 
detector elements for detecting the rays of the 
beam, the at least one focus is rotated relative to 
the object under examination on at least one focal 

10 path running around the object under examination 

with opposite detector array, the detector elements 
of the detector array supplying initial data which 
represents the attenuation of the rays as they pass 
through the object under examination, 

15 - the initial data is filtered, 

- the filtered initial data is then backpro j ected 
three-dimensionally in order to produce at least one 
slice of a layer of the object under examination 
which has a layer thickness, the slice representing 

20 absorption values, obtained from the initial data, 

of the voxels belonging to the layer of the object 
under examination for the radiation of the beam, and 

- during the backpro j ection, the rays are weighted as a 
function of their position in the beam. 

25 

In a particular embodiment, the filtering is carried 
out in the direction of the tangent to the focal path 
belonging to the respective focal position. It has been 
shown that, with this filtering direction, a 

3 0 particularly high image quality can be achieved. The 
selection of this filtering direction is based on the 
finding that the abovedescribed 2D method permitting a 
high image quality and based on preliminary images 
would, so to speak, change into a 3D method if the 

3 5 sections of the focal path on which the calculation of 
preliminary images is based were to be shortened so 
severely that they would cover only a single projection 
whose data would then be filtered in the direction of 
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the tangent to the focal path, and in that it would 
then be possible to expect that such a 3D method would 
then permit a comparatively good image quality as the 
2D method. 

5 

If the beam has an extent in the direction of rotation 
and an extent in the direction of the axis of rotation, 
then it is advantageous to weight rays which are 
arranged centrally in the beam, as based on the extent 
10 of the beam in the direction of the axis of rotation, 
to a greater extent than the rays which are arranged 
close to the edge in the beam, as based on the extent 
of the beam in the direction of the axis of rotation. 

15 The method according to the invention may then be 
carried out particularly easily if/ before filtering, a 
conversion of the initial data obtained in fan beam 
geometry in the form of rays P(a,(3,q) into parallel 
data present in parallel beam geometry in the form of 

20 rays P(6,f3,q) (azimuthal "rebinning") or P(9,p,q) 
(complete "rebinning" , that is to say azimuthal and 
radial "rebinning") is performed. In this case, with 
reference to fig. 3 



25 a is the focal angle 

(3 is the fan angle 

qr is the row index of the detector system 

corresponding to the z coordinate, 

0 = a + (3 is the parallel fan angle, 

30 p = R F sin({3) is the parallel coordinate corresponding 

to the distance of the ray from the axis 
of rotation (system axis), and 

R F is the radius of the focal path. 



3 5 According to a preferred embodiment of the invention, 
the backprojection of the parallel data is carried out 
in that, in the course of the backprojection for each 
voxel V(x,y, z) , for each 6e [o,^] for the rays 
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P(9 + kn,p,q) and PiO + kn.p.q) whose projection along the 
system axis goes through (x,y), the sum 



is formed, where 

x,y,z are the coordinates of the respective voxel 
V(x,y, z) , 

k is a whole number corresponding to the number 

10 of half revolutions of the focus included in 

the reconstruction, 
p are the parallel coordinates of those rays 

whose projections along the system axis run 
through the coordinates (x,y) of the respective 
15 voxel V(x,y,z), 

P are the fan angles of those rays whose 

projections along the system axis run through 
the coordinates (x,y) of the respective voxel 
V(x, y, z) , 

2 0 h is a weighting function determining the layer 

thickness of the layer of the object under 
examination represented in the slice produced, 
d is a function which is equal to the distance of 

the respective ray from the corresponding voxel 
25 V(x,y) and or is dependent on the distance of 

the respective ray from the corresponding voxel 
V(x,y) , and 

W represents a weighting function which weights 

rays with a large parallel fan angle 9 less 

3 0 than rays with a small parallel fan angle 6. 

In this case, the notation j~j expresses the fact that 

the formation of the sum can optionally be carried out 
for rays obtained by means of azimuthal "rebinning" or 
35 complete "rebinning", the filtering tangential to the 
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focal path in the case of azimuthal "rebinning" being 
filtering in the j8 direction and, in the case of 
complete ^rebinning", being filtering in the p 
direction . 

5 

Because of the summing both over k and over q, it is 
then ensured that all the rays running through one and 
the same voxel are taken into account and the radiation 
dose supplied to the object under examination is thus 
10 used completely. 

A particularly preferred embodiment provides that, for 
the purpose of the backpro j ection of the parallel data, 
the sum 

15 




normalized to the sum H of the weights h 



20 



k q 



0 + kn 9 



25 



is formed. This procedure permits an image quality 
which is improved again, since possible overemphasis of 
voxels which are struck by more rays than other voxels 
is eliminated, and thus corresponding artefacts are 
avoided. The CT value of the respective voxel is 
obtained by summing over 0. 



According to the invention, the weighting function can 
3 0 be a function of the parallel fan angle with W(d+kn) , 
and preferably represent a smooth function which has 
the value 1 for central rays, as based on the extent of 
the fan of rays in the direction of the axis of 
rotation or z direction, and tends to 0 for marginal 
3 5 rays. 
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If the detector array has detector elements arranged in 
the manner of rows, then the weighting can be 
represented by a function of the row number W(q) , a 
5 smooth function also preferably being. used here which, 
for rays with a centrally placed detector row or a 
plurality of centrally placed detector rows, has a high 
value, preferably the value 1, and for rays arranged at 
the margin tends to 0; for example this can be a cos 2 
10 function. 

The method according to the invention can be .used, for 
example when, according to one variant of the 
invention, the focal path is a circular path (tomogram 

15 scanning) . According to a preferred variant of the 
invention, however, the focal path is a spiral path, 
which is brought about by the focus being moved on a 
circular path about the system axis and, at the same 
time, there being a relative movement between focus and 

20 object under examination in the direction of the system 
axis. On the basis of such spiral scanning, even 
relatively large volumes of the object under 
examination can be examined without problems. 

2 5 In the case of tomogram scanning, k is normally k=l or 

k=2 ; in the case of spiral scanning, k is selected 
while taking account of the relative displacement in 
the direction of the system axis carried out per full 
revolution, such that the region to be imaged of the 

3 0 object under examination is registered completely. 

According to the basic idea of the invention, the 
inventors also propose to improve a CT device for 
scanning an object under examination with a beam 
35 originating from at least one focus and with a flat 
detector array having a large number of distributed 
detector elements for detecting the rays of the beam, 
the at- least one focus moving relative to the object 
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under examination on at least one focal path running 
around the object under examination with the opposite 
detector array, to the effect that at least means for 
collecting detector data, filtering and backpro j ection 
5 are provided, which carry out the method described 
above, these means preferably and as far as possible 
being implemented by programs or program modules . 

The invention will be explained in more detail below 
10 using an exemplary embodiment which is illustrated in 
the appended schematic drawings, in which: 

Fig. 1 shows a CT device having a plurality of rows of 
detector elements in a partly perspective, 
15 partly block-diagram illustration; 

Fig. 2 shows a longitudinal section through the device 
according to Fig. 1; 

Fig. 3 shows a diagram illustrating "rebinning"; 

Fig. 4 shows scanning with opposite focus and detector 
20 in the beam with a rotating gantry; 

Fig. 5 shows the course of two different weighting 
functions W; 

Fig. 6 shows scanning with opposite focus and detector 
in the beam with the gantry rotating and moving 
25 in the z direction. 



Figs 1 and 2 illustrate a CT device of the third 
generation suitable for carrying out the method 
according to the invention. Its measurement 

3 0 arrangement, designated overall by 1, has an x-ray 
source, designated overall by 2, with a beam aperture 
stop 3 arranged in front of the latter, close to the 
source (Fig. 2) and a detector system 5 formed as a 
two-dimensional array of a plurality of. rows and 

35 columns of detector elements - one of these is 
designated 4 in Fig. 1 - with a beam aperture stop 6 in 
front of the latter, close to the detector (Fig. 2). In 
Fig. 1, for clarity, only eight rows of detector 
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elements 4 are shown. However, as indicated by dots in 
Fig. 2, the detector system 5 has further rows of 
detector elements 4. 

5 The x-ray source 2 with the aperture stop 3, on the one 
hand, and the detector system 5 with the aperture stop 
6, on the other hand, are fitted, in the manner which 
can be seen from Fig. 2, to a rotary frame 7, opposite 
each other in such a way that a pyramidal x-ray beam 

10 originating from the x-ray source 2 during operation of 
the CT device and masked by the adjustable aperture 
stop 3, whose marginal rays are designated 8, strikes 
the detector system 5. In this case, the aperture stop 
6 is adjusted in accordance with the cross section of 

15 the x-ray beam set by the aperture stop 3 such that 
only that region of the detector system 5 which can be 
struck directly by the x-ray beam is- exposed. In the 
operating mode illustrated in Figs 1 and 2, these are 
eight rows of detector elements 4, which are referred 

20 to as active rows in the following text. The further 
rows, indicated by dots, are covered by the aperture 
stop 6 and therefore not active. 

Each row of detector elements 4 has a number K of 
2 5 detector elements, (3 k = fix to (3 K being the channel index 
and each detector element being assigned a fan angle 
(3 k . The fan angle of the central detector element is 
equal to 0; the fan angles of the two outermost 
detector elements are /3i = +/3 /nax and (3 K = - /3 max . 

30 

The active rows L Q of detector elements 4 are 
designated L 2 to L Q , q = 1 to Q being the row index 
which, in the case of the exemplary embodiment 
described, corresponds so to speak to the z coordinate. 

35 

The x-ray beam has the cone angle <p plotted in Figs 1 
and 2, which is the opening angle of the x-ray beam in 
a plane containing the system axis Z and the focus F. 
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The opening angle of the x-ray beam in a plane at right 
angles to the system axis Z and containing the focus F 
(fan opening angle) is 2/3 JT7ax and is plotted in Fig. 1. 

5 The rotary frame 7 can be set rotating by means of a 
drive device 22 about a system axis designated by Z. 
The system axis z runs parallel to the z axis of a 
three-dimensional rectangular coordinate system 
illustrated in Fig. 1. 

10 

The columns of the detector system 5 likewise run in 
the direction of the z axis, while the rows, whose 
width b is measured in the direction of the z axis and 
is 1 mm, for example, run transversely with respect to 
15 the system axis Z and the z axis. 

In order to be able to move an object under 
examination, for example a patient, into the beam path 
of the x-ray beam, a mounting device 9 is provided, 

20 which can be displaced parallel to the system axis Z, 
that is to say in the direction of the z axis, to be 
specific in such a way that there is synchronization 
between the rotational movement of the rotary frame 7 
and the translational movement of the mounting device, 

25 with the effect that the ratio between translation and 
rotational speed is constant, this ratio being 
adjustable by a desired value for the advance v of the 
mounting device per revolution of the rotary frame 
being selected . 

30 

It is therefore possible for a volume of an object 
under examination located on the mounting device 9 to 
be examined in the course of a volume scan, it being 
possible for the volume scan to be carried out in the 
3 5 form of spiral scanning, with the effect that, with 
simultaneous rotation of the measuring unit 1 and 
translation of the mounting device 9, a large number of 
projections from different projection directions is 
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recorded by means of the measuring unit for each 
revolution of the measuring unit 1. During the spiral 
scanning, the focus F of the x-ray source moves 
relative to the mounting device 9 on a spiral path 
5 designated S in Fig. 1. The spiral scanning must 
extend in the a direction over at least 7t+2(i maXf in order 
to permit the complete reconstruction of a CT image for 
each row of detector elements, but can also be longer 
as desired within the technical limits of the CT 
10 device. 

However, because of the fact that a plurality of rows 
of detector elements 4 are present, a volume of the 
object under examination can also be examined in the 

15 course of what is known as a tomogram scan, in which no 
relative movement takes place in the direction of the z 
axis between the measuring unit 1 and mounting device 9 
(v = 0). In the case of the tomogram scan, the size of 
the volume examined is therefore determined by the 

20 number of active rows of detector elements 4. During a 
tomogram scan, the focus F moves on a circular focal 
path which lies in a plane designated the mid-plane 
below. 

2 5 The tomogram scanning can be carried out in the form of 

a part revolution or in the form of a complete 
revolution, the part revolution covering a part 
revolution interval of at least 7r+2p max (half a 
revolution plus the fan opening angle) , which permits 
30 complete reconstruction of a CT image, while a full 
revolution covers 2k. 

The measured data read out in parallel from the 
detector elements of each active row of the detector 

3 5 system 5 during the spiral or tomogram scanning and 

corresponding to the individual projections P(cx,p>,q) in 
fan beam geometry is subjected to digital /analog 



200210627 - 12 - 

conversion in a data conditioning unit 10, serialized 
and transmitted to an image computer 11. 



Following preconditioning of the measured data in a 
5 preconditioning unit 12 belonging to the image computer 
11, the resultant data stream passes to a slice 
reconstruction unit 13, which uses the- measured data to 
reconstruct slices of desired layers of the object 
under examination in accordance with a method according 
10 to the invention based on "filtered backpro j ection" and 
yet to be described in detail. 

The CT images are composed of pixels - (pixel = picture 
element) assembled in the manner of a matrix, the 

15 pixels being assigned to the respective image plane, 
each pixel being assigned a CT number in Houns field 
units (HU) and the individual pixels being displayed in 
accordance with a CT number/gray value scale in a gray 
value corresponding to their respective CT number. In 

20 this case, each pixel illustrates a voxel (voxel 

volume element) of the layer illustrated in the CT 
image and belonging to the object under examination. 
Since, because of the multiplicity of rows of the 
detector system 5 and, if appropriate, the spiral 

25 scanning, measured data with respect to a plurality of 
layers of the object under examination is obtained, 3D 
data is available, which is subjected to 3D 
backpro j ection within the context of the invention. As 
the final result, 3D image data is available in the 

30 form of a three-dimensional matrix, for example with 
the axes x,y, z, each element in the matrix 
corresponding to a voxel V(x,y,z) and containing the 
gray value corresponding to the associated CT number. 
Those elements of the three-dimensional matrix which 

3 5 have the same ^c, y, or z value then in each case 
represent a planar slice of the layer of the object 
under examination corresponding to the defining x, y, 
or z value. 
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The images reconstructed by the slice reconstruction 
unit 13 are displayed on a display unit 16, for example 
a monitor, connected to the image computer 11. 

5 

The x-ray source 2, for example an x-ray tube, is 
supplied with the necessary voltages and currents, for 
example the tube voltage U, by a generator unit 17. In 
order to be able to set these to the respectively 
10 necessary values, the generator unit i7 is assigned a 
control unit 18 with keyboard 19, which permits the 
necessary adjustments . 

In addition, the other operation and control of the CT 
15 device is carried out by means of the control unit 18 
and keyboard 19, which is illustrated by the fact that 
the control unit 18 is connected to the image computer 
11 . 

20 Inter alia, the number Q of active rows of detector 
elements 4, and therefore the position of the aperture 
stops 3 and 6, can be adjusted, for which purpose the 
control unit 18 is connected to adjusting units 2 0 and 
21 assigned to the aperture stops 3 and 6 . 

2 5 Furthermore, the rotation time r needed by the rotary 

frame 7 for a complete revolution can be adjusted, 
which is illustrated by the fact that the drive unit 22 
assigned to the rotary frame 7 is connected to the 
control unit 18 . 

30 

Although it is possible in principle to implement the 
method according to the invention in fan beam geometry 
as well, the CT device described is preferably operated 
in a mode in which the method according to the 

3 5 invention is implemented in parallel beam geometry. 

Accordingly, the data obtained during the scanning of 
the region of the body of the patient which is relevant 
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to the respective examination by means of spiral or 
tomogram scanning in fan beam geometry is first of all 
converted, in a manner known per se, into data in 
parallel beam geometry by means of a method generally 
5 designated "rebinning" . This conversion is based on 
resorting the data obtained in fan beam geometry in 
such a way that rays are removed from different 
projections recorded in fan beam geometry, and joined 
together to form a pro j ection in parallel beam 
10 geometry. In parallel beam geometry, data from an 
interval with a length n is sufficient to be able to 
reconstruct a complete image. In order to obtain this 
data, data in fan beam geometry from an interval of 
length n+2p max must nevertheless be available. 

15 

A projection in parallel beam geometry is illustrated 
in Fig. 3. Accordingly, all n parallel rays RPi to RP N 
of this projection assume the parallel fan angle 6 with 
respect to the x axis of the coordinate system which is 

2 0 illustrated in Fig. 3 and coincides with that according 

to Fig. 1. 

By using the parallel ray RPi illustrated by a 
continuous line in Fig. 3, the transition from fan beam 
25 to parallel been geometry will be explained below. 

The parallel ray RPi originates from a pro j ection 
obtained in fan beam geometry for the focal position Fi 
located on the focal path S . The central ray RF Z i 

3 0 belonging to this projection in fan beam geometry, 

running through the axis of rotation 14 and therefore 
the z axis of the coordinate system, is likewise 
plotted in Fig. 3. The focal position Fi corresponds to 
the focal angle Ofi; this is the angle formed by the x 
35 axis and the central ray RF zi . As compared with the 
central ray RF z i, the ray RPi has the fan angle /3. It is 
thus easy to see that, for the parallel fan angle 9 it 
is true that 6 = a + (3 . 
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The beam spacing p from the axis of rotation 14 or the 
z axis, measured at right angles to the respective 
parallel ray, is given by p=R F sin((5) . 

5 

As becomes clear by using the central ray RP Z 
illustrated in a thickened line in Fig. 3 and extending 
through the axis of rotation 14 and the x axis, this 
ray is the central ray of a projection in fan beam 
10 geometry recorded in fan beam geometry for the focal 
position F z at the focal angle a z - Since it is true that 
/?=0 for the central ray of a projection recorded in fan 
beam geometry, it becomes clear that the following is 
true of the case of central rays: 

15 

Depending on whether an azimuthal or complete 
"rebinning" is carried out, the parallel projections 
are present in the form 
P(a, 

20 or in the form 

ro,p, q) 

where 

a is the focal angle 

(3 is the fan angle 

25 q is the row index of the detection system 

corresponding to the z coordinate, 
9=a+{3 is the parallel fan angle 

p=R F sin ( (3) is the parallel coordinate corresponding to 
the ray spacing from the axis of rotation 
30 (system axis), and 

R F is the radius of the focal path. 



In an operating mode corresponding to a first 
embodiment of the method according to the invention, 
35 which can be selected by means of the keyboard 19, the 
CT device described operates on the basis of 
projections obtained by azimuthal "rebinning" . The data 
corresponding to these projections is filtered in the (3 
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direction, that is to say in each case in the direction 
of the tangent T belonging to the focal position of the 
central ray of the respective parallel projection (see 
fig. 3), specifically by using one of the filter 
5 kernels customary in computer tomography, e.g. a Shepp- 
Logan or Ramachandran-Lakshminarayanan kernel . 

The parallel data filtered in this way is then 
backprojected in that, in the course of the 
10 backproj ection, for each voxel V(x,y,z) for each #e[0,7r] 
for the rays P(& + k7r,/3,q) whose projection along the 
system axis goes through (x,y), the sum 



15 

is formed, where 

x,y, z are the coordinates of the respective voxel 
V(x,y, z) , 



k q 



20 



k 



is a whole number corresponding to the number 
of half revolutions of the focus included in 
the reconstruction, 



p 



are the fan angle as of those rays whose 
projections along the system axis run through 
the coordinates (x,y) of the respective voxel 
V(x, y, z) , and 



25 



h 



is a weighting function which determines the 

layer thickness of the layer of the object 

under examination illustrated in the slice 
produced, and 



30 



is a function which is equal to the distance of 
the respective ray from the corresponding voxel 
V(x,y) or is dependent on the distance of the 
respective ray from the corresponding voxel 
V(x,y), and 



35 



W(q) 



represents a weighting function which weights 
rays which are positioned peripherally or at 
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least decentrally in the beam, as based on the 
extent in the direction of the axis of 
rotation, less than rays which are arranged 
centrally. 

5 

As a result of the selected filtering direction and as 
a result of the summation both over k and over q, 
firstly cone beam artefacts are avoided, in the 
interest of a high dose utilization, all the rays 
10 running through a voxel V(x,y,z) are taken into 
account, but in addition, the positioning of the ray in 
the beam is also taken into account by means of 
appropriate weighting. 

15 The absorption value Vx,y,z associated with a voxel 
V(x,y,z) is obtained by means of summation over 6 over 
at least half a revolution, that is to say by forming 

M x,y,z j ^x,y,z (^0 

20 

The CT number corresponding to the respective 
absorption value is determined in a conventional way 
from the absorption value. 

25 In this case, different weighting functions h and 
different functions d can be set by means of the 
keyboard 19 . 

For example, a triangular or trapezoidal function is 
3 0 suitable as a weighting function h. 

The weighting function W is represented as a function 
of the projection angle 6. 

3 5 The function d can be the distance of the respective 
parallel ray from the voxel V(x,y,z) or, instead, for 
example the z (axial) component of this distance. 
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In a modification of the first operating mode described 
above, for the purpose of backpro j ection of the 
parallel data, the sum 

5 

0 

normalized to the sum H of the weights h 

k q 

10 

is formed. This permits an image quality which is again 
improved, since the possible overemphasis of voxels 
which are "illuminated", that is struck by rays, in a 

15 plurality of half revolutions is eliminated and thus 
corresponding artefacts are avoided. This redundancy 
occurs during spiral scans when the relative 
displacement taking place per full revolution of the 
measuring arrangement is so low (low pitch) that voxels 

2 0 are irradiated repeatedly. 

A second operating mode, which corresponds to a further 
embodiment of the method according to the invention and 
can be selected by means of the keyboard 19, differs 

2 5 from the first operating mode in that the CT device 
described does not operate on the basis of projections 
obtained by means of azimuthal but by means of complete 
"rebinning" . The data corresponding to these 
projections is filtered in the p direction in the case 

30 of complete "rebinning", that is to say in each case 
likewise in the direction of the tangent T belonging to 
the focal position of the central ray of the respective 
parallel projection (see Fig. 3). 
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Accordingly, for the parallel data filtered in this way 
in the course of the backpro j ection, the sum 



P,y, ( 0 ) = 77 2 X Kb) ■ h i d *.y, fa + A 4 P fa + A </) 

5 

is formed, /5 being the parallel coordinates of those 
rays whose projections along the system axis run 
through the coordinates (x,y) of the respective voxel 
V(x,y, z) . 

10 

In the case of the second operating mode, too, in a 
modification relating to the backpro j ection of the 
parallel data, a sum normalized to the sum H of the 
weights h 

15 

p , y , (*) = Z Z w ia) ' h (<>, (* + k *> p> 9% p{o + P, q) 

k q 

is formed, namely the sum 

ft q 

20 

In the case of the first and second operating modes 
described above, a function of the method according to 
the invention is provided in which, in conjunction with 
25 a voxel V(x,y,z), all the rays whose projection along 
the axis of rotation of 14 or the z axis goes through 
x,y are considered. Whether and to what extent these 
rays are taken into account is determined by the 
weighting functions W and h and the function d. 
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However, the CT device can also have further operating 
modes which can be selected by means of the keyboard 19 
and which correspond to those described previously, 
with the difference that, for a given focal position, 
5 the theoretical ray running through the respective 
voxel V(x,y,z) is determined and then, taking into 
account the weighting function h and the function d in 
the sum formation, only those rays which can actually 
make a contribution to the sum are included in the 
10 course of the backpro j ection . 

In the • exemplary embodiments described, the relative 
movement between the measuring unit 1 and the mounting 
device 9 is in each case produced by the mounting 

15 device 9 being displaced. However, within the scope of 
the invention, there is also the possibility of 
arranging for the mounting position 9 to be fixed in 
position and, instead, to displace the measuring unit 
1. In addition, within the scope of the invention, 

20 there is the possibility of producing the necessary 
relative movement by displacing both the measuring unit 
1 and the mounting device 9. 

In Figure 4, the problem of weighting the rays in a 2D 
25 representation is shown, simplified ■ and illustrated 
schematically. In this example, the gantry of the CT 
with focus and detector in the 0° position is shown by 
continuous lines and in the 180° position, rotated 
about the z axis by 180° , is shown by dots. The 
30 detector has 12 detector rows Ll to L12, which are each 
struck by the associated rays SI to S12 from the beam. 
It should be pointed out that the spatial dimensions 
shown are illustrated exaggerated with respect to the 
extent of the detector, in order to be able to explain 
35 the invention more clearly. The gantry shown moves on 
a circular path here about the z axis or system axis. 
Movement in the direction of the z axis does not take 
place . 
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If the voxel V2 which lies centrally in the beam is 
considered, this voxel will be penetrated by the 
central rays S4 in both positions of the gantry. 
5 According to the invention, the measurements from these 
rays aire highly weighted. At the same time, in the 0° 
position, the voxel VI which lies somewhat more 
peripherally is penetrated in the 0° position by the 
outer ray SI and by the ray S4 lying further in. In 

10 accordance with the weighting according to the 
invention, the ray SI will be given less weight than 
the ray S4 . If the likewise peripherally arranged voxel 
V3 is considered, then the weighting here is reversed. 
Since, here, the detector is constructed with detector 

15 elements which are divided up into rows and columns in 
the manner of a matrix, the weighting of the rays can 
be coupled directly to the row number q. 

In figure 5, located below, three examples of a course 
2 0 according to the invention of weightings W q are 
illustrated. Common to all is the greater weighting of 
the central rays with respect to the outer rays. 
However, the functions are different. The dotted line 
represents a course which drops off linearly on either 
25 side of the centre, the dashed line results in 
identical weighting over a large central region with a 
tendency to drop off rapidly at the periphery, and the 
continuous line shows a bell-like course which permits 
a gentle transition between the extreme values 0 and 1. 

30 

Figure 6 shows CT beams from a gantry of a spiral CT 
respectively in the 0° position and adjacent 180° 
position. In accordance with the essence of the spiral 
CT, these positions are offset in relation to each 
35 other in addition to the rotation in the z direction. 
In this illustration, the voxel VI is penetrated by the 
ray SI both in the 0° position and in the 18 0° 
position. In accordance with this decentral position 
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of the ray SI and the associated detector row Ll 
located on the outside in the detector, these two 
measurements are weighted lowly. By contrast, the 
voxel V2 in the 0° position is penetrated by the ray 
5 S12 located decentrally, which is weighted lowly in 
this position for the voxel V2 , but in the 180° 
position is penetrated by the central ray S6, which is 
weighted highly for this voxel V2 . In accordance with 
the position of the rays SI in the 0° and 180° 
10 positions, these rays are both weighted lowly in 
relation to the voxel VI. Overall, as a result of this 
additional weighting W illustrated here, the result is 
a substantial improvement in the image quality. 

It should be pointed out that the conical x-ray beam 
has a rectangular cross section in the exemplary 
embodiment described here, but other cross-sectional 
geometries are also possible within the scope of the 
invention . 

In connection with the exemplary embodiments described 
above, . CT devices of the third generation are used, 
that is to say the x-ray source and the detector system 
are displaced jointly about the system axis during the 
image generation. However, the invention can also be 
used in connection with CT devices of the fourth 
generation, in which only the x-ray source is displaced 
about the system axis and interacts with a stationary 
detector ring, if the detector system is a multi-row 
array of detector elements. 

The method according to the invention can also be used 
in CT devices of the fifth generation, that is to say 
CT devices in which the X-radiation originates from a 
35 plurality of foci of one or more x-ray sources 
displaced about the system axis, and not just from one 
focus, if the detector system has a multi-row array of 
detector elements . 



20 
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The CT devices used in connection with the exemplary 
embodiments described above have a detector system of 
detector elements arranged in the manner of an 
5 orthogonal matrix. However, the invention can also be 
used in connection with CT devices whose detector 
system has a flat array with detector elements arranged 
in another way. 

10 The exemplary embodiments described above relate to the 
medical application of the method according to the 
invention. However, the invention can also be applied 
outside medicine, for example in checking luggage or in 
material examination . 
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Patent claims 



1. A method of creating images in computer 
tomography, comprising the method steps: 

5 1.1 to scan an object under examination with a beam 
originating from at least one focus and with a 
flat detector array having a large number of 
distributed detector elements for detecting the 
rays of the beam, the at least one focus is 

10 rotated relative to the object under examination 

on at least one focal path running around the 
object under examination with opposite detector 
array, the detector elements of the detector array 
supplying initial data which represents the 

15 attenuation of the rays as they pass through the 

object under examination, 

1.2 the initial data is filtered, 

1.3 the filtered initial data is then backpro j ected 
three-dimensionally in order to produce at least 

20 one slice of a layer of the object under 

examination which has a layer thickness, the slice 
representing absorption values, obtained from the 
initial data, of the voxels belonging to the layer 
of the object under examination for the radiation 

2 5 of the beam, and 

1.4 during the backpro ject ion, the rays are weighted 
as a function of their position in the beam. 

2. The method as claimed in claim 1, characterized in 
30 that the filtering is carried out in the direction 

of the tangent to the focal path belonging to the 
respective focal position. 

3. The method as claimed in either of claims 1 and 2, 

3 5 characterized in that the beam has an extent in 

the direction of rotation and an extent in the 
direction of the axis of rotation, rays which are 
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arranged centrally in the beam, as based on the 
extent of the beam in the direction of the axis of 
rotation, being weighted to a greater extent than 
the rays which are arranged close to the edge in 
5 the beam, as based on the extent of the beam in 

the direction of the axis of rotation. 



4. The method as claimed in one of claims 1 to 3 , 
characterized in that, before filtering, a 

10 conversion of the initial data obtained in fan 

beam geometry in the form of rays P(a,(i,q;) into 
parallel data present in parallel beam geometry in 
the form of rays P(9,(3,q) or P(0,p,q) is 
performed, where 

15 a is the focal angle 

(3 is the fan angle 

gr is the row index of the detector 

system corresponding to the z 
coordinate, 

20 6 = a + (3 is the parallel fan angle, 

p = R F sin({3) is the parallel coordinate 

corresponding to the distance of the 
ray from the axis of rotation (system 
axis ) , and 

2 5 R F is the radius of the focal path. 

5. The method as claimed in claim 4, characterized in 
that the backpro jection of the parallel data is 
carried out in that, in the course of the 

30 backpro jection for each voxel V(x,y,z), for each 6 

g [o,7i] for the rays P(0 + kn,(3,q) and P(6 + kn,p,q) 

whose projection along the system axis goes 
through (x, y) , the sum 



• w W=ZZ^-^( tf+ **.{|}.«))- / { tf+ * ,r -{|}'9 



35 



200210627 - 26 - 

is formed, where 

x,y,zare the coordinates of the respective voxel 
V(x,y, z) , 

k is a whole number corresponding to the 

5 number of half revolutions of the focus 

included in the reconstruction, 
p are the parallel coordinates of those rays 

whose projections along the system axis run 
through the coordinates (x,y) of the 
10 respective voxel V(x,y,z), 

J3 are the fan angles of those rays whose 

projections along the system axis run 
through the coordinates (x,y) of the 
respective voxel V(x,y,z), 
15 h is a weighting function determining the 

layer thickness of the layer of the object 
under examination represented in the slice 
produced, 

d is a function which is equal to the distance 

2 0 of the respective ray from the corresponding 

voxel V(x,y) or is dependent on the distance 
of the respective ray from the corresponding 
voxel V(x,y) , and 
W represents a weighting function which 

25 weights rays with a large parallel fan angle 

9 less than rays with a small parallel fan 
angle 9. 

6. The method as claimed in claim 5, characterized in 

3 0 that, during the backpro j ection of the parallel 

data, the sum 




35 



normalized to the sum H of the weights . h 
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is formed. 

5 7. The method as claimed in either of claims 5 and 6, 
characterized in that the weighting function 
represents a function of the parallel fan angle 
with W(0+Iar). 

10 8. The method as claimed in claim 7, characterized in 
that the weighting function W(Q+krf) represents a 
smooth function which has the value 1 for 
centrally arranged rays in the beam, based on the 
parallel fan angle, and tends to 0 for rays 

15 arranged at the edge. 

9. The method as claimed in either of claims 4 and 5, 
characterized in that the detector array has 
detector elements arranged in the manner of rows, 

2 0 and the weighting function represents a function 

of the row number W(q) . 

10. The method as claimed in claim 8, characterized in 
that the weighting function W(g) represents a 

2 5 smooth function which has the value 1 for rays to 

at least one centrally located detector row and 
tends to 0 for rays to detector rows at the edge. 

11. The method as claimed in one of claims 1 to 9, 
30 characterized in that the focal path is a circular 

path. 

12. The method as claimed in one of claims 1 to 9 , 
characterized in that the focal path is a spiral 

3 5 path which is brought about by the focus being 

moved about the system axis on a circular path 
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and, at the same time, a relative movement between 
focus and object under examination in the 
direction of the system axis taking place. 

5 13. The method as claimed in one of claims 1 to 12, 
characterized in that the detector elements on the 
detector array are arranged distributed in rows 
and lines, preferably in the manner of a matrix. 

A CT device for scanning an object under 
examination, with a beam originating from at least 
one focus and with a flat detector array having a 
large number of distributed detector elements for 
detecting the rays of the beam, the at least one 
focus moving relative to the object under 
examination on at least one focal path running 
around the object under examination with the 
opposite detector array, at least means for 
collecting detector data, filtering and 
backpro jection being provided as claimed in one of 
claims 1 to 13 . 



10 14. 



15 



15 The CT device as claimed in claim 14, 
characterized in that the aforementioned 
25 functional means are at least partly implemented 

by programs or program modules . 
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Abstract 

Method of creating images in computer tomography, and 
CT device for carrying out the method 

The invention relates to a method for image 
reconstruction for computer tomography with a flat, 
extended detector, the rays being weighted during the 
backpro j ection as a function of their position in the 
beam . 

In addition, the invention relates to a CT device for 
carrying out the above-mentioned method. 
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